o 
O 



> 



oo 



>$ 



Mid-Infrared Spectra of Be Stars 

S. A. Rinehart} J. R. Houck, & J. D. Smith 
Cornell University 

rinehart@tristan.tn. cornell.edu, jrhl3@cornell.edu, 
jdsmithOastrosun . tn . Cornell . edu 



0\ ■ ABSTRACT 



We present the first medium-resolution (R ~ 600) mid-infrared (8-13. 3/xm) spectra 
of 11 Be stars. A large number of lines are observed and identified in these spectra, 
including, as an example, 39 hydrogen recombination lines in the spectrum of 7 Cas. 
In the majority of our spectra, all of the observed lines are attributable to hydrogen 
recombination. Two of the sources, /3 Lyr and MWC 349 also show emission from other 
10 i species. Both of these objects show evidence of [Ne II] emission, and /? Lyr also shows 

evidence of He I emission. We tabulate the effective line strength and line widths for the 
observed lines, and briefly discuss the physical implications of the observed line series. 
We also use a simple model of free-free emission to characterize the disks around these 



0> ■ sources. 

Mh| Subject headings: Infrared Radiation — Stars: Emission-line, Be 

o 



1. Introduction 



Be stars are defined as B stars which display optical emission in one or more Balmer lines of 
F^ ', hydrogen, or which have displayed such emission at some point in the past. There are also several 

subdivisions of Be stars, including the "classical" Be stars and Herbig Be stars. Herbig Ae/Be stars 
are considered the intermediate mass (5-20 Mq) counterparts to T Tauri stars; they are pre- main 
sequence stars surrounded by an optically thick, dusty disk remaining from proto-stellar collapse. 
From this disk arises forbidden line emission and Ha emission. Classical Be stars, on the other 
hand, are intermediate mass main-sequence objects, believed to have an extended optically thin 
gaseous envelope. While these two types of objects share similar nomenclature, they have little in 
common outside of their masses. 

The sources presented in this paper are all classical Be stars, and henceforth the term Be star 
refers to classical Be stars. Classical Be stars are also divided into two categories, "non-shell". 
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or "normal" stars, and "shell" stars. Shell stars display broad emission wings in one or more of 
the hydrogen lines and narrow absorption bands for ionized metals (Slettebak 1988), features not 
observed in non-shell Be stars. Several suggestions about the relation between shell stars and normal 
Be stars have been made: this relation is complicated by the fact that several stars, including 7 Cas 
and Pleione (HR 1180), have been observed to change from "normal" spectra to "shell" spectra 
(Doazan & Thomas 1982). Hanuschik (1996) explains shell stars as Be stars which are observed 
equator-on, with the transition between normal and shell spectra explained by density variations 
in a concave equatorial disk viewed at a small inclination angle (~15°). Alternatively, it has been 
suggested that the shell phase is a common evolutionary phase of Be stars (Slettebak, et al. 1992), 
involving the formation of a gas shell outside the stellar envelope. The formation of the gas shell 
can explain the transition from normal to shell spectra, while the dissipation of the gas shell causes 
the transition from shell to normal spectra. 

Optical spectroscopy has raised a number of questions about Be stars. Individual spectral lines 
are found to be double peaked, with a ratio of the violet peak flux to the red peak flux (the V/R 
ratio) which differs from unity. Further, this ratio and the line strengths of the hydrogen Balmer 
lines have been found to be time- variable in many Be stars (Telting et al 1993). 

Near-infrared (NIR) spectroscopy has added to this puzzle. Variability of line strengths and 
V/R ratios of the double-peaked lines mirror observations in the optical (Hamann & Simon 1989). In 
addition, a number of the hydrogen recombination lines which appear in the NIR have anomalously 
high line strengths, and the NIR continuum has been observed to vary with time (Dougherty & 
Taylor 1994). 

Ultraviolet observations complicate the problem further. The strong lines observed in the UV 
should arise from a fast, low-density wind (Smith 1995). This is at odds with the optical and NIR 
observations, which require a relatively dense gas to produce the observed features. A possible 
explanation for this discrepancy is discussed below. 

X-ray observations have shown several Be stars to be X-ray binaries. 7 Cas, the brightest of 
the northern hemisphere Be stars has a compact companion, most likely a white dwarf (Haberl 
1995). Other Be stars are also in binary systems; a large number of Be stars, including many of 
the sources presented in this sample, are spectroscopic binaries, including /3 Lyr, C, Tau, and (j) Per 
(Jared et al. 1989). 

Of all of the wavelength regions, the mid-infrared has been the least explored. Photometry 
by Gehrz et al (1974) from 2.3 to 19.5/im showed that the spectral energy distribution of the IR 
excess was consistent with free-free emission from a warm disk-shaped envelope. It has also been 
suggested (when discussing NIR observations) that studies of infrared line emission can be very 
useful in probing the stellar envelope of Be stars, as the emission is highly dependent upon the 
temperature and density structure of the envelope (Marlborough et al. 1997). This is one of the 
motivations for the study introduced here. 

Explaining the multitude of observations by modelling is a daunting task. The most common 



model used to explain the observations is the axisymmetric equatorial disk model of Poeckert & 
Marlborough (1978). A slow, high-density wind at the equator creates a thick envelope which 
produces the observed optical and NIR emission, while the UV emission arises from a fast, tenuous 
wind originating at the poles of the star (Lamers & Waters 1987). The presence of an equatorial 
disk has been supported by interferometric and polarimetric observations (Quirrenbach et al. 1997, 
Stee et al. 1995), but these models are incomplete, as they still cannot explain some Be star 
phenomenae, such as the observed V/R variability. 

A more complete discussion of Be stars can be found in several excellent review articles, 
including Slettebak (1988) and Dachs (1986). 

In this paper, we present the first high-sensitivity, medium resolution {R ^ 600) mid-infrared 
(8-13. 3^m) spectra of be stars. In §2, we discuss the observations and data reduction techniques. 
In §3, we present the observed spectra, and briefly mention characteristics of particular sources. 
In §4, we discuss using free-free emission from the stellar envelope to model the MIR continuum. 
Finally, in §5, we briefly discuss some of the implications of these observations. 



2. Observations and Data Reduction 

The observations presented here were taken using SCORE (van Cleve et al. 1998) at the 
Cassegrain focus of the Hale 5-m telescope.^ SCORE is a novel medium resolution mid-infrared 
spectrograph built originally as a proof of concept of the Infrared Spectrograph (IRS) short-hi 
module for SIRTF. It uses cross-dispersion to obtain a 7.5-15/um spectrum in a single exposure. 
The obtained spectrum has a resolution of ~600, with a 1" x 2" slit. This short slit length is 
ideal for sources which are unresolved in the MIR at Palomar. SCORE also makes use of a second 
MIR array for slit-viewing, which both simplifies source acquisition during observing and provides 
information for photometric calibration (see discussion below) . A listing of all of the observations is 
presented in Table 1. All of these observations were made using standard beam-switched chopping 
and nodding techniques. Typical integration times ranged from 240 to 1090 seconds, and the entire 
N-band was observed at a single setting for each object. The la sensitivity at ll.O^m in 100 
seconds of integration for our instrument is roughly 90mJy. The MIR seeing for these observations 
was typically 1". 

The individual integrations of each object were stacked, and the N-band spectra were extracted 
using SCOREX, extraction software designed to optimally extract the spectrum from the cross- 
dispersed format of SCORE (Smith et al. 1998). These spectra were then calibrated using observed 
spectra of standard stars (Cohen, et al. 1992) to properly remove spectral features from the sky 
background. The observed spectra from these standard stars were taken at nearly the same airmass 
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as our observed objects, as even relatively small differences in airmass result in different strengths of 
atmospheric features, leading to inaccurate subtraction of background features from our calibrated 
spectra. 

Absolute photometry for the spectra was found by using the slit-viewer images of SCORE. 
Since a small chop amplitude (~5") was used for these observations, the chopped images appear in 
the field-of-view of the slit-viewer. We calculate the slit-throughput by comparing the flux incident 
on the slit-viewer focal plane with the total flux found in the observed spectrum; mathematically, 

_ Fi,(Slitviewer) 

~ f F,{\)g{X) dX ^> 

where A is the renormalization factor, Fu{X) is the flux as a function of wavelength on the spectro- 
graph, g{X) is the slit viewer filter passband function, and F,^(Slitviewer) is the total flux incident on 
the slit-view focal plane array. The renormalization factor A should account for the slit throughput. 
To provide the absolute calibration, we use the same procedure for our calibrator star, deriving an- 
other normalization factor, Acaiib- By multiplying the observed, calibrated spectrum of our source 

by -J we are able to normalize the observed spectrum to have the correct amount of flux, without 

introducing errors from instrumental effects. The absolute normalization of the spectra is uncertain 
to roughly 5%, due primarily to the uncertainty in the slit throughput efficiency. 



3. Observed Spectra of Be stars 

In Figure 1, we present the reduced spectra of our objects. These spectra clearly show the 
presence of a large number of spectral lines. Further, in most of the spectra, there are a number 
of lines which are weak relative to the strong MIR continuum. Since these weak lines appear in 
several different spectra, but not in calibration star spectra, we conclude that they are not artifacts. 
Further, the clear presence of multiple lines in each series implies the presence of other lines within 
the series. For instance, the H 29^10 is weak in the spectrum of 7 Cas, but the presence of other 
transitions, (including the H 29^9 transition in this case) support the identification of this weak 
line, since they arise from a common excited state. Similar support can be found for nearly all 
of the weak observed lines. Finally, because all of the lines occur at wavelengths attributable to 
hydrogen recombination (for most of our objects, see below), there is strong evidence to support 
the pure-hydrogen recombination spectra. 

There is also another test for the presence of these weak lines. Nine of the 11 spectra display 
similar strong hydrogen line emission (all but (3 Lyr and MWC 349), but a number of these spectra 
have low signal to noise. These nine sources will hereafter be referred to as the hydrogen spectra Be 
stars (or HS stars). If we suppose that all of the lines observed in the high signal to noise spectra 
(such as 7 Cas and C, Tau) are present in the spectra of all HS sources, we can coadd the spectra 
from all nine of these objects. If these lines are, in fact, present in a majority of the sources, we 



expect the weak lines to be accentuated. If, on the other hand, these weak hnes are present in only 
a few spectra, we expect them to be diluted (relative to their appearance in the individual spectra). 
This coadded spectrum is shown in Figure 2, with the line identifications marked at the top of the 
plot. Examination of this spectrum shows a large number of lines, including a number which are 
only observed as very faint peaks in the spectra of our bright sources. This supports the idea that 
all of the observed hydrogen transitions are present in the majority of the HS sources. 

In Table 2, we list all of the observed lines, their line widths and effective line strengths, 
and the hydrogen recombination transition responsible for their production, in order of decreasing 
wavelength. This table includes the data for all of the HS sources (which excludes /? Lyr and 
MWC 349). The data from the two peculiar objects {/3 Lyr and MWC 349) is presented in Table 3, 
because of the clear differences between their spectra and the spectra of the majority of our sources. 
Note that in Table 2, every observed line is attributable to hydrogen recombination. We examined 
the spectra for the presence of line species other than hydrogen, in particular for the presence of 
He I or He II, since these species can produce emission lines coincident with hydrogen lines. He 
I is easily identified by the presence of the 10.88 /xm He I 3S-3Po transition, while He II can be 
identified by the presence of many additional hydrogen-like transitions. In the nine HS spectra, 
there is no evidence for the He I 10.88//m transition or for the numerous He II transitions which 
would be present if the species accounted for any of the observed emssion. 

In P Lyr, however, we observed the He I 3S-3Po transition at 10.88 /um. Since Helium re- 
combinations occur at the same wavelength as hydrogen recombinations (for instance, both H 9— >7 
and He I 9— >7 occur at 11.3//m) the presence of this line implies that the observed lines are not 
necessarily purely due to hydrogen, as is true for most of our other sources. The possible overlap 
of hydrogen and helium recombination line makes exact identification of these lines problematic. 
We also find evidence for the [Ne II] 12.8)Um transition, a line also observed in the spectrum of 
MWC 349. MWC 349 showed no evidence for He I or He II emission, and only showed a few weak 
hydrogen transitions (with the exception of the H 7— >6 transition which was quite strong). We will 
return to a discussion of the observed spectral features of these objects later in this section. 



3.1. Hydrogen Spectra Be Stars 

While the nine HS sources presented here have similar MIR spectra, they are significantly 
different types of objects in several ways. Seven of the sources are classified as shell stars. Five 
of our sources are spectroscopic binaries, and three of the sources are X-ray binaries. Further, 
while most Be stars do not display radio emission, four of the sources presented here are relatively 
bright radio sources. The characteristics of the individual objects are tabulated in Table 4. For 
completeness, some of the characteristics of the two peculiar sources are also listed on this table. 

The large differences in the properties of the individual objects leads us to look for correlations. 
Are any of these individual properties related to the MIR observations presented here? We looked 



at correlation plots for each of the different categories listed in Table 4. An example of one of 
these correlation plots is shown in Figure 3. We have plotted the ratio of the H 9^7 to H 7^6 
transition versus the ratio of the H 10^7 to H 7^6 transitions in this figure, with the XRB plotted 
as + symbols, spectroscopic binaries as open triangles, and single stars as open diamonds. We 
have examined a number of similar correlation plots, looking for evidence that our MIR spectra are 
effected by each of these properties, but have found no such evidence. However, given the small 
statistics (only nine sources), this lack of detection is not surprising, and only by building a more 
significant set of data will such correlation tests be valuable for understanding the effects of these 
peculiarities on our MIR spectra. 

The presence of so many high-level hydrogen transitions provides valuable insight into the origin 
of line emission. The line strengths are inconsistent with optically thin line emission (Hummer & 
Storey 1987), and therefore must originate at optical depths of ~1. The optically thick emission will 
simply be the product of the Planck function, the line width, and the surface area of the emission 
region. Since the Planck function of a gas with temperatures consistent with the production of 
such high H recombination series decreases with wavelength and our observations show no clear 
trend in line width with wavelength, the surface area must be increasing with wavelength. ISO 
observations have shown that the line width actually decreases with wavelength (Hony et al. 1999), 
strengthening this statement. This implies a density gradient, in order to balance line ratios; if the 
density is uniform, then the low-level transitions will become optically thick much more quickly 
than their high-level counterparts, giving them an effective emission surface significantly larger 
than the high-level lines, and overcompensating for the decreasing strength with wavelength. This 
was previously noted from NIR observations of 7 Cas (Hamann &: Simon 1987). These authors 
concluded that gradients as small as r~^ were sufficient to explain their observed line ratios. 



3.2. Peculiar Be Stars 

P Lyr and MWC 349 display spectra which are significantly different from our HS sources, and 
from each other. Each of these sources is unique, and it is therefore not surprising that their MIR 
spectra should be peculiar as well. We briefly discuss the peculiar aspects of the MIR spectra of 
these objects here, and present some possible explanations for these peculiarities. 

P Lyr is a known spectroscopic and eclipsing binary. The two members of this binary are B 
stars, with masses of Mi ~ 13Mq and M2 ~ 3Mq (Harmanec & Scholz 1993). The two stars are in 
close orbit, with the mass-loss from the less-massive star filling its Roche lobe, overflowing to form 
an accretion disk around the more massive star (De Greve & Linnell 1994). This situation is quite 
different from the disk-like stellar envelopes around our HS Be stars, which occupy a comparatively 
small volume. 

The accretion process can also possibly explain observed He I emission by invoking convective 
dredging in the donor star, thus bringing helium into the Roche lobe and the accretion disk. The 



presence of the He recombination lines in this spectrum, when compared to the HS sources is 
most easly explained by either a higher He abundance in the circumstellar disk or a higher level of 
excitation. If we compare (3 Lyr to 7 Cas, however, we see that 7 Cas is significantly hotter than 
P Lyr {Tpiyr < 20000K). Since no He emission is observed in 7 Cas, whose disk gas is at higher 
excitation, we conclude that the observed He emission in /? Lyr is due to a greater He abundance. 

The biggest question raised from examination of the MIR spectrum of (3 Lyr is the cause of 
[Ne II] emission. While a great number of emission lines have been observed in the optical and 
NIR spectra of /3 Lyr, there has been no evidence for forbidden emission in any of these spectra 
(Johnson 1978), or for mid-infrared emission from forbidden lines in our spectra ([Ar III] at 9.0;um, 
[S IV] at 10.5/um). Why should [Ne II] emission, an no other forbidden emission, be observed? 
The two properties which should be examined in particular are the ratio of the observed to critical 
density (ne/ncrit) for [Ne II] emission and the energy required to form Ne II. Upon examination, 
we find that both of these characteristics of Ne are favorable for the observed forbidden emission. 
[Ne II] has a very high critical density (10^ cm~^), higher than any other astronomically strong 
MIR forbidden transition. Further, it has a relatively low ionization potential (21.5 eV), making it 
possible to ionize a large enough fraction of the Neon gas by the hot B stars to create the observed 
emission. 

Using the observed strength of the [Ne II] line, the distance to /3 Lyr (270 pc) from the 
Hipparchos catalog (1997), and making a few assumptions about the emitting gas (e.g. Ug = 
Ucrit, Tf, = lO^K, solar abundances), we can calculate an approximate mass of emitting gas. This 
calculation gives a mass of the emitting gas of Mtot ~ 10^^ Mq. This compares favorably with 
calculations of the mass-transfer rate from the "donor" star to the "accretor" star, if we assume 
that some small fraction of mass-transfer escapes into an extended gas shell. Such calculations, 
assuming complete mass conservation between the two stars, give values of M ~ 2 x 10^^ M© yr^" 
(Hubeny et al. 1994). If we assume that the gas is fully ionized, we can estimate an emission volume 
for the [Ne II] line. We assume, for this approximation, that the emission volume is spherical, and 
find that the emitting region should be of order 100 AU across. For comparison, the separation 
of the two components of the binary is only of order 20 AU. From these estimates, we suggest 
that the forbidden line emission does not arise in the Roche lobe or accretion disk, but must come 
from a much more extended gas shell around the binary. This could also explain the difference 
in the observed hydrogen spectra of this source, relative to our HS sources. The extended, low 
density outer gas shell will produce optically thin line emission from hydrogen, while the disk 
produces optically thick line emission. The contribution from the optically thin emission will alter 
the hydrogen line ratios, explaining why the line ratios for hydrogen emission from f3 Lyr are not 
~1 (as for our HS sources). 

MWC 349 has the largest MIR to optical flux ratio of all of our sources, with this ratio more 
than 100 times larger than our other sources. This could be due to several different phenomena, 
including a much cooler source or much more circumstellar mass (including, perhaps, a large amount 
of dust, leading to severe reddening). Independent of the cause, the red color of MWC 349, 



compared to our HS sources, implies that the nature of this object or its environment may be 
quite different from classical Be stars. Observations of MWC 349 in the submillimeter detected Ha 
transitions which were greatly amplified; MWC 349 was the first source observed with hydrogen 
laser emission (Martin-Pintado et al. 1989). Further, recent ISO observations have detected more 
than 160 emission lines between 2.4 and 190//m, including every hydrogen a transition (H5 -^ 4 
to H16^ 15) in this range (Thum et al. 1998), showing the presence of infrared lasers from the 
large amplification of the Ha lines. Further, lasing/masing has been observed from a number of 
the hydrogen /3 lines. It is thought that the conditions which allow the observed masing/lasing 
from MWC 349 include the high temperature of (and corresponding large ultraviolet flux from) the 
central star, a dense, massive, neutral disk (much more massive than the disks in most Be stars), 
and the coincidental edge-on view we have of the disk (Strelnitski et al. 1996). 

A significant amount of forbidden line emission in the optical has been observed in MWC 349 
(Brugel & Wallerstein 1979, Andrillat & Swings 1976, Allen & Swings 1976)). We have also detected 
forbidden emission, in the form of the [Ne II] 12.81/um line. It is unlikely that this forbidden 
emission arises from the same region as the hydrogen emission, however. Recent work by Them & 
Greve (1997) used the observed Paschen decrement to estimate the electron density in the disk at 
He = 10^, over one hundred times greater than the critical density for [Ne II] forbidden emission. 
The forbidden emission more likely arises in an extended gas component around the star. Such an 
extended component has been observed via radio measurements and has been associated with a 
slow (50 km s^^) gaseous outflow. Cohen et al. (1985), using observations from the VLA, found 
that the radio observations were well-matched by a spherical 1/r^ wind model, with mass-loss rates 
of order IO^^Mq yr^^ and a temperature of roughly 9000 K. 

Using the same technique as described for (3 Lyr, we estimate the mass of the emitting gas 
from the forbidden emission. We use the same assumptions as in the case of (3 Lyr, and estimate 
a distance to the source of 400 pc (Thompson et al. 1977). From these, we find a mass of the 
emitting gas of Mtot ~ 1O~^M0. The mass of the disk should be substantially larger than this, up 
to several solar masses (Thompson et al., Thum &: Greve 1997). However, the forbidden emission 
we observe will only come from the regions around the star where the density is relatively low, and 
therefore will only arise from a small fraction of the disk, since the majority of the disk has a high 
(~ lO^cm-3) density (Thum et al. 1998). 



4. Continuum Emission 

The majority of photospheric radiation from hot stars is at short wavelengths, in the optical 
or ultraviolet (an 16000 K star has a peak flux at ~ 1800A, well into the UV). In the MIR, 
photospheric emission is well inside the Raleigh-Jeans tail of the blackbody distribution, producing 
relatively weak emission in this wavelength range. Classical Be stars have been long known to 
produce large IR excesses, and a common explanation for this is emission via free-free processes. 



Free-free radiation from ionized gas is a powerful source of emission at long wavelengths. 
Typically, such emission has been used to explain radio emission from cool astronomical sources 
(Panagia & Felli 1975), but it can also explain the infrared excesses observed from many sources. Be 
stars have well-known IR excess, and free-free emission from the extended stellar envelope provides 
a convenient explanation for this excess. 

Gehrz et al. (1974) used models of optically thin and optically thick free- free emission to match 
infrared photometry of Be stars from 2.3//m to 19.5/im. Waters et al (1984) successfully modeled 
IRAS Be star observations with an optically thick free-free emission. We have used models of free- 
free emission to fit our MIR continuum observations, finding values for the shell radius {Rsh)-, the 
electron density (ng), and the shell temperature (Tg/j). 



4.1. Calculations of Free- Free Emission 

Three emission processes are included in our model of the MIR continuum: emission from the 
photosphere of the star, optically thin free-free emission, and optically thick free-free emission. 

Because neither of the free- free processes are very effective at short wavelengths, we can con- 
strain the stellar parameters without having to consider the effects of the shell. In the infrared, the 
photospheric spectrum is essentially a blackbody produced at the photospheric temperature, so the 
Planck function describes this emission. Initial values for the temperature (T^) and radii (i?^) of 
our sources were approximated from their respective spectral types (Waters et al. 1987), and we 
then adjusted the values of both T^ and i?^ to agree with J-band (1.25^m) and K-band (2.20^m) 
photometry (Gezari et al. 1993). Distances to the stars were obtained from the Hipparchos Catalog 
(1997). All of these parameters, including the J-band and K-band fiuxes, are listed in Table 5. 

To fit the observed MIR continuum, we require several assumptions. First, we assume that the 
MIR emission arises from an extended stellar envelope around the star. Further, it is assumed that 
this envelope is flattened into an oblate spheroidal disk, with a semi-minor axis ~ 1/10-Rsh where 
Rsh is the semi-major axis of the spheroid. Further, we assume that the stellar envelope has both 
uniform density and uniform temperature. These assumptions are clearly not correct, but greatly 
simplify calculations without introducing significant amounts of error. Finally, we assume that we 
are looking into the disk-like envelope edge-on. This last assumption is supported by two facts: 
the large values of v sin i for these sources, and the central line absorption of shell spectra (seven 
of our nine HS sources are shell stars, see Table 4). 

The MIR emission from the stellar photosphere is calculated using the Planck function (making 
use of i?*, T^, and D), as described above, but is corrected by applying an extinction factor in the 
form of e~'^^P^ "^' . This will account for free-free absorption of photospheric radiation as a function 
of wavelength. Thus, we write 
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Next, we include the effect of optically thick free-free emission. The optically thick emission 
will be characterized by the Planck function at the shell temperature (Tgh)- This emission will arise 
from an area of the disk which is defined by the surface a depth p into the disk, where the optical 
depth is 
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This cuts off the optically thick emission at short wavelengths, as the absorption coefficient becomes 
very small and the entire disk becomes optically thin. Mathematically, 
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where x{t) is the ratio of the area of the optically thick emission to the maximum area of optically 
thick emission. 

Finally, we include the effect of the optically thin region of the emission disk. This can be 
accounted for to first order by integrating the emission function for free- free emission {el ) over the 
volume of the stellar envelope. We improve upon this first-order calculation by also including the 
extinction (e^'^^''^) factor in the integral over volume. In this way, we account for the reduction in 
optically thin emission with increasing wavelength. This gives us, replacing Equation 3, 



To fit our observations, we sum the emission from these three sources, varying the temperature 
(Tsh), radius (Rsh), and electron density (ue) of the shell. This sum is then compared to our 
observed spectrum. We calculate a reduced x^ value based upon nine points in our observed 
spectra; these nine points located in regions at least O.OS^m from any of the observed emission 
lines, in order to not contaminate the x^ with line emission. The flux at each point is determined 
by averaging the flux over A — 0.05 to A + 0.05^m, in order to reduce the effects of local noise spikes. 
The results of these fits are shown in Table 6. In Figure 4, we show the observed spectrum of 7 Cas 
with the overlaid fit. The observed spectrum is shown as a dotted line, the fit curve is shown as a 
solid line, the dashed line and dot-dash line are the optically thin and optically thick components 
of the emission respectively, the diamonds are the nine points used to calculate the reduced x^^ 
and the open triangle is the 12;um flux measured by IRAS. 
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4.2. Comparisons with Previous Work 

Two works in particular, Gehrz et al. 1974 and Waters et al. 1987, have calculated parameters 
of Be stellar envelopes. Gehrz et al. calculated Rgh and Ue from their observations, based upon an 
estimate of Tg}-^ (this estimate was constrained solely by the requirement that T^^ be between lOOOOK 
and T^) by using the estimated Tgh and by calculating the wavelength where the extrapolations 
of optically thin and optically thick free-free emission intersect (at which point the optical depth 
r(Ac) = 1). This provides a relation between the radius of the disk and the electron density. 
Further, by assuming a Rayleigh-Jeans distribution for both the optically thick shell flux and the 
stellar flux (at Ac), they obtain an expression for the radius of the shell, which then determines the 
electron density of the disk. For these calculations, they have assumed the disk to be flattened with 
an aspect ratio A = Tg/ddisk = 5 (where ddisk is the full thickness of the disk), the same value used 
for our modelling. As in our model, theirs includes the simplifying assumption that the electron 
density was uniform in the disk. 

Waters et al. assume Tgh = 0.8T^ and calculate Rgh and the density profile {p{r)) for the 
stellar disks using the "curve of growth" method (Lamers & Waters 1984). For this method, they 
assume a disk with an opening angle 6, a disk density which is proporational to r^", a density at 
the inner edge of the shell of po, and a disk radius r^/j (the sharp outer cutoff to the disk). For 
their modelling, they have assumed a pole-on view of the disk, such that the opening angle 6 is 
deprecated in their calculations; this is contrary to the both our assumption and the assumption 
of Gehrz et al. that the disk is viewed edge-on. They then use a combination of optically thin and 
optically thick free-free radiation as a comparision to the observed monochromatic infrared excess 
(as found from IRAS spectra and using assumed stellar parameters). By comparing the shapes of 
the observed excess and their models, they obtain a value fo the density parameter n, and from the 
shift required to match the observed and model curves they obtain po. These models are dependent 
upon the assumed shell radius, such that the authors run a series of models to find the best-fit 
combination of n and r^/i, yielding the values used for comparision in this paper (notably, their 
values for Vgh are all lower limits). We do not compare our values for n^ with the p{r) calculated by 
Waters et al. because of the large uncertainty in the mass densities they derived, the uncertainty in 
the proper ionization fraction of the envelope gas, and the dependence upon the exact shell radius. 

We present the results of the modelling discussed in this paper, as well as the results of both 
Gehrz et al. and Waters et al. in Table 7 for ease of comparision. We find that our results are 
generally consistent (within factors of a few) with those of both authors. However, in a number of 
the cases. Waters et al. derive shell radii which are significantly larger than the radii derived here. 
This can be explained by the assumptions which have gone into our modelling, and the wavelength 
regimes of the observations. We have assumed an isothermal disk, whereas Waters et al. did not. 
Since their observations utilize IRAS data, they will be sensitive to the cool parts of the shell which 
should exist in the case of a thermal gradient. These discrepancies are particularly apparent in the 
cases of (j) Per and 7 Gas, bright sources with accurate 12;um, 25^m, and 60^m IRAS detections. 
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Comparing our results to those of Gehrz et al, we find that in the case of C Tau we predict an 
electron density which is a factor of two below that of Gehrz et al, and a moderately larger disk. 
In the cases of 7 Cas and (p Per, we derive larger, hotter disks than Gehrz et al, but have similar 
electron densities. The simplicity of the models used by Gehrz et al. are most likely responsible for 
this discrepancy. They extrapolated curves for both the optically-thin and optically-thick segments 
of the free-free emission, then used the intercept point of these two curves to derive numerical values 
for radius and temperature of the shell. The improved modelling techniques used here should be 
better able to limit these same parameters. 

In Table 7, we also show calculated total mass of hydrogen in the emitting disk. These 
numbers assume complete ionization of the hydrogen gas, such that rig = nn, and uniform density, 
an assumption already used in the modelling. We find that the masses of the disks are relatively 
smaU, ranging from 2.5xlO~-'^^M0 (o Aqr) to 8.1xlO~^M0 (7 Cas). 



5. Discussion 

These spectra present MIR features of Be stars and give some indications of the range of 
emission characteristics of them as well. Continuum emission is explained by a combination of 
optically-thick and optically-thin free-free emission. Both the continuum and line emission is likely 
to arise in a warm stellar envelope. The peculiar objects {(3 Lyr and MWC 349) display complicated 
line emission which is difficult to explain with simple models. In the case of (3 Lyr, a relatively 
simple explanation for the presence of hydrogen and helium is possible, and the observed [Ne II] 
emission is easily explained as emission from a large gas cloud surrounding the binary system. In 
the case of MWC 349, the observed features are consistent with previous observations, and the 
observed forbidden emission implies a large volume of low-density ionized gas. 

The authors would like to thank the staff at Palomar Observatory for their assistance. We 
would also like to thank S. Hony for sharing ISO data and results prior to publication. SAR 
acknowledges L.B.F.M. Waters and M. Simon for helpful discussions. This research was partly 
supported by NASA Contract 960803. 
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Table 1: Observations of Be Stars 



Object 


Date 


Integration 


Spectral 


IRAS 






Observed 


Time (s) 


Type 


12;um Flux 


(Jy) 


V' Per 


Aug. 29, 1998 


668 


B5Ve 






Tj Tau 


Aug. 29, 1998 


653 


B7III 


4.32 




C Tau 


Aug. 28, 1998 


481 


B2IV 


8.30 




48 Per 


Aug. 28, 1998 


661 


B3Ve 


3.80 




(/>Per 


Aug. 27, 1998 


741 


B2Vpe 






7 Cas 


Aug. 27, 1998 


240 


BOIVe 


18.83 




Aqr 


Aug. 27, 1998 


586 


B7IVe 


1.18 




/3Lyr 


July 5, 1998 
July 13, 1998 


930 


B7Ve+ 


5.01 




EWLac 


July 12, 1998 


661 


B3IVpe 






K Dra 


July 5, 1998 


1094 


B6IIIpe 


4.12 




MWC 349 


July 12, 1998 
Aug. 27, 1998 


481 


Bpe 


179.0 
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Table 2A. Data on observed lines - ip Per 



Transition Xtheory Kbs AA Fj, 



H 14- 


-^9 


12.587 


12.582 ± 0.0064 


13.40 ± 5.833 


6.73 ± 3.96 


H 19- 


-^10 


12.611 


12.613 ± 0.0024 


4.06 ± 2.622 


3.04 ± 2.93 


H 11- 


-^8 


12.387 


12.391 ± 0.0047 


3.55 ± 4.034 


2.81 ± 4.21 


H7^ 


6 


12.372 


12.371 ± 0.0016 


7.18 ± 1.236 


10.60 ± 2.43 


H21- 


-^10 


11.792 


11.787 ± 0.0046 


8.93 ± 4.136 


4.14 ± 2.56 


H 15- 


-^9 


11.540 


11.540 ± 0.0042 


13.10 ± 4.098 


7.83 ± 3.09 


H9^ 


7 


11.309 


11.308 ± 0.0013 


82.50 ± 1.216 


9.86 ± 1.91 


H24- 


-^10 


11.033 


11.033 ± 0.0049 


6.91 ± 4.302 


2.82 ± 2.24 


H25- 


-^10 


10.855 


10.854 ± 0.0041 


5.78 ± 3.545 


2.55 ± 1.89 


H 16- 


-^9 


10.804 


10.803 ± 0.0018 


7.52 ± 1.640 


7.09 ± 2.03 


H26- 


-^10 


10.701 


10.695 ± 0.0053 


8.00 ± 1.414 


3.03 ± 1.67 


H 12- 


-^8 


10.504 


10.502 ± 0.0021 


8.11 ± 1.920 


5.97 ± 1.90 


H 17- 


-^9 


10.261 


10.256 ± 0.0032 


12.90 ± 2.921 


7.87 ± 2.40 


H30- 


-^10 


10.258 




Blended with 17- 


-^9 


H 18- 


-^9 


9.847 


9.842 ± 0.0035 


8.42 ± 3.154 


6.40 ± 3.09 


H 19- 


-^9 


9.522 


9.518 ± 0.0043 


6.66 ± 3.782 


3.31 ± 2.59 


H 13- 


-^8 


9.392 


9.390 ± 0.0013 


6.11 ± 1.192 


9.35 ± 2.29 


H20- 


-^9 


9.261 


9.258 ± 0.0051 


8.63 ± 4.691 


4.39 ± 3.10 


H21- 


-^9 


9.047 


9.044 ± 0.0046 


4.35 ± 3.810 


2.36 ± 2.57 


H 10- 


-^7 


8.760 


8.758 ± 0.0012 


4.97 ± 1.373 


7.77 ± 2.78 


H 14- 


-^8 


8.665 


8.663 ± 0.0022 


8.41 ± 1.997 


10.34 ± 3.23 


H24- 


-^9 


8.594 


8.594 ± 0.0033 


4.61 ± 3.466 


3.75 ± 3.47 
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Table 2B. Data on > 


observed lines - r] Tau 




Transition 


■^theory 


■^obs 


AA 


F, 






(urn.) 


(fim) 


(10--Vm) (10-- 


^9 W cm-2 ^m-i) 


H 14- 


-^9 


12.611 


12.586 ± 0.0032 


9.14 ± 3.075 


5.07 ± 2.2.8 


H 11- 


-^8 


12.387 


12.393 ± 0.0016 


4.21 ± 1.782 


2.80 ± 1.74 


H7^ 


6 


12.372 


12.375 ± 0.0010 


5.68 ± 1.392 


6.76 ± 2.07 


H 15- 


-^9 


11.540 


11.538 ± 0.0026 


8.84 ± 2.376 


5.27 ± 1.86 


H9^ 


7 


11.309 


11.308 ± 0.0009 


7.41 ± 0.870 


8.67 ± 1.33 


H23- 


-^10 


11.243 


11.245 ± 0.0070 






H24- 


-^10 


11.033 


11.044 ± 0.0093 






H25- 


-^10 


10.855 


10.859 ± 0.0042 


6.02 ± 5.527 


1.75 ± 2.17 


H 16- 


-^9 


10.804 


10.802 ± 0.0018 


7.74 ± 1.621 


6.50 ± 1.79 


H 12- 


-^8 


10.504 


10.501 ± 0.0016 


7.65 ± 1.483 


5.61 ± 1.43 


H 17- 


-^9 


10.261 


10.257 ± 0.0028 


10.00 ± 2.586 


4.97 ± 1.73 


H 18- 


-^9 


9.847 


9.844 ± 0.0026 


6.72 ± 2.346 


4.43 ± 2.05 


H 19- 


-^9 


9.522 


9.516 ± 0.0026 


5.23 ± 2.019 


3.10 ± 1.65 


H 13- 


-^8 


9.392 


9.389 ± 0.0019 


3.14 ± 2.605 


2.62 ± 2.89 


H 10- 


-^7 


8.760 


8.758 ± 0.0016 


6.42 ± 1.429 


9.55 ± 2.80 


H 14- 


-^8 


8.665 


8.663 ± 0.0023 


7.19 ± 2.113 


8.19 ± 3.22 


H24- 


-^9 


8.594 


8.595 ± 0.0084 






H 15- 


-^8 


8.155 


8.164 ± 0.0017 


5.57 ± 1.506 


7.33 ± 2.65 


H29- 


-^9 


8.173 




Blended with 15^8 
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Table 2C. Data on ( 


Dbserved lines - C Tau 




Transition 


■^theory 


^obs 


AA 


F, 






(^m) 


(fim) 


(10-Vm) (10-: 


^9 W cm-2 ^m-i) 


H 18- 


-^10 


13.188 


13.191 ± 0.0025 


7.09 ± 2.144 


4.41 ± 1.76 


H 19- 


-^10 


12.611 


12.611 ± 0.0036 


8.54 ± 3.164 


5.16 ± 2.27 


H 14- 


-^9 


12.587 


12.586 ± 0.0032 


6.88 ± 3.412 


4.23 ± 2.42 


H 11- 


-^8 


12.387 


12.395 ± 0.0062 


3.80 ± 9.182 


0.91 ± 2.40 


H7^ 


6 


12.372 


12.375 ± 0.0039 


11.60 ± 3.078 


8.93 ± 2.89 


H20- 


-^10 


12.157 


12.160 ± 0.0027 


7.94 ± 2.588 


4.52 ± 1.94 


H21- 


-^10 


11.792 


11.796 ± 0.0030 


10.40 ± 2.794 


5.05 ± 1.76 


H 15- 


-^9 


11.540 


11.542 ± 0.0025 


10.00 ± 2.219 


6.20 ± 1.83 


H22- 


-^10 


11.492 


11.493 ± 0.0040 


8.41 ± 3.562 


2.93 ± 1.67 


H9^ 


7 


11.309 


11.310 ± 0.0022 


8.51 ± 2.032 


5.46 ± 1.73 


H23- 


-^10 


11.243 


11.243 ± 0.0032 


7.54 ± 2.775 


3.53 ± 1.72 


H24- 


-^10 


11.033 


11.037 ± 0.0037 


6.13 ± 4.111 


2.56 ± 2.31 


H25- 


-^10 


10.855 


10.860 ± 0.0037 


7.93 ± 2.941 


3.95 ± 1.93 


H 16- 


-^9 


10.804 


10.804 ± 0.0016 


9.71 ± 1.481 


10.54 ± 2.12 


H26- 


-^10 


10.701 


10.696 ± 0.0036 


4.66 ± 4.525 


2.08 ± 2.53 


H27- 


-^10 


10.567 


10.573 ± 0.0058 


4.18 ± 5.796 


0.88 ± 1.57 


H 12- 


-^8 


10.504 


10.507 ± 0.0039 


6.73 ± 3.551 


3.74 ± 2.41 


H28- 


-^10 


10.451 


10.458 ± 0.0080 


6.84 ± 7.271 


1.91 ± 2.48 


H29- 


-^10 


10.348 


10.353 ± 0.0113 






H 17- 


-^9 


10.261 


10.264 ± 0.0016 


7.96 ± 1.505 


8.71 ± 2.18 


H30- 


-^10 


10.258 




Blended with 17^9 




H31- 


-^10 


10.177 


10.178 ± 0.0024 


3.84 ± 2.404 


1.89 ± 1.65 


H 18- 


-^9 


9.847 


9.849 ± 0.0016 


7.23 ± 1.451 


9.15 ± 2.39 


H 19- 


-^9 


9.522 


9.527 ± 0.0010 






H 13- 


-^8 


9.392 


9.395 ± 0.0019 


7.21 ± 1.708 


9.64 ± 3.01 


H20- 


-^9 


9.261 


9.266 ± 0.0018 


5.88 ± 1.719 


7.34 ± 2.84 


H21- 


-^9 


9.047 


9.049 ± 0.0036 


6.44 ± 3.288 


5.61 ± 3.65 


H22- 


-^9 


8.870 


8.875 ± 0.0035 


5.78 ± 2.963 


4.66 ± 3.22 


H 10- 


-^7 


8.760 


8.761 ± 0.0024 


5.72 ± 1.965 


5.45 ± 3.03 


H23- 


-^9 


8.721 


8.723 ± 0.0025 


6.64 ± 2.232 


7.86 ± 3.46 


H 14- 


-^8 


8.665 


8.667 ± 0.0024 


6.96 ± 2.207 


8.30 ± 3.48 


H24- 


-^9 


8.594 


8.596 ± 0.0031 


5.84 ± 2.774 


6.67 ± 4.02 


H25- 


-^9 


8.485 


8.488 ± 0.0021 


2.75 ± 2.296 


3.02 ± 3.51 


H26- 


-^9 


8.391 


8.392 ± 0.0033 


4.99 ± 2.427 


4.60 ± 3.26 


H27- 


-^9 


8.309 


8.306 ± 0.0036 


4.37 ± 3.477 


3.55 ± 3.58 


H28- 


-^9 


8.236 


8.239 ± 0.0011 






H 15- 


-^8 


8.155 


8.156 ± 0.0009 


4.85 ± 1.035 


11.02 ± 3.38 


H29- 


-^9 


8.173 




Blended with 15^8 




H30- 


-^9 


8.155 


8.103 ± 0.0023 


6.33 ± 2.102 


8.12 ± 3.74 
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Table 2D. Data on observed lines - 48 Per 



Transition 


■^theory 


Ao6s 


AA 


F, 






(^m) 


(^m) 


(10-3//m) (10 


-19 w cm-2 ^m-i) 


H 18- 


-^10 


13.188 


13.188 ± 0.0065 


6.14 ± 4.530 


3.04 ± 2.95 


H 19- 


-^10 


12.611 




Blended with H 14- 


^9 


H 14- 


-^9 


12.587 


12.587 ± 0.0028 


9.00 ± 2.685 


6.27 ± 2.42 


H 11- 


-^8 


12.387 




Blended with H 7^ 


■6 


H7^ 


6 


12.372 


12.379 ± 0.0019 


8.33 ± 1.594 


7.72 ± 2.01 


H20- 


-^10 


12.157 


12.159 ± 0.0033 


5.90 ± 3.344 


2.30 ± 1.66 


H21- 


-^10 


11.792 


11.792 ± 0.0061 


10.10 ± 5.602 


2.66 ± 1.92 


H 15- 


-^9 


11.540 


11.537 ± 0.0027 


8.52 ± 2.482 


4.18 ± 1.60 


H9^ 


7 


11.309 


11.306 ± 0.0011 


6.41 ± 0.889 


6.00 ± 1.12 


H23- 


-^10 


11.243 


11.236 ± 0.0062 


7.72 ± 5.487 


1.69 ± 1.60 


H24- 


-^10 


11.033 


11.028 ± 0.0044 


4.16 ± 8.134 


0.95 ± 2.39 


H 16- 


-^9 


10.804 


10.802 ± 0.0019 


7.31 ± 1.777 


4.82 ± 1.56 


H27- 


-^10 


10.567 


10.566 ± 0.0197 






H 12- 


-^8 


10.504 


10.501 ± 0.0015 


6.79 ± 1.334 


5.05 ± 1.35 


H28- 


-^10 


10.451 


10.446 ± 0.0226 






H 17- 


-^9 


10.261 


10.257 ± 0.0013 


7.61 ± 1.179 


6.69 ± 1.37 


H 18- 


-^9 


9.847 


9.846 ± 0.0021 


6.19 ± 1.812 


4.92 ± 1.94 


H 19- 


-^9 


9.522 


9.517 ± 0.0019 


4.36 ± 1.390 


3.65 ± 1.54 


H 13- 


-^8 


9.392 


9.388 ± 0.0011 


5.92 ± 0.913 


8.08 ± 1.64 


H20- 


-^9 


9.261 


9.260 ± 0.0049 


9.45 ± 4.608 


4.07 ± 2.60 


H21- 


-^9 


9.047 


9.044 ± 0.0046 


6.15 ± 4.097 


3.14 ± 2.68 


H22- 


-^9 


8.870 


8.870 ± 0.0051 


2.58 ± 5.069 


1.33 ± 2.77 


H 10- 


-^7 


8.760 


8.757 ± 0.0009 


5.82 ± 0.756 


9.98 ± 1.70 


H23- 


-^9 


8.721 


8.718 ± 0.0059 


5.76 ± 5.464 


2.14 ± 2.60 


H 14- 


-^8 


8.665 


8.662 ± 0.0015 


5.68 ± 1.286 


7.55 ± 2.27 



20 







Table 2E. Data on ^ 


observed lines - 


(/>Per 




Transition 


'^theory 


^ohs 


AA 




Fu 






(/wm) 


(^m) 


(10~3/im) 


(10- 


19 w cm-2 ^m-i) 


H 18- 


-^10 


13.188 


13.185 ± 0.0059 


13.10 ± 5.435 




4.23 ± 2.29 


H 19- 


-^10 


12.611 


12.615 ± 0.0009 


3.84 ± 1.099 




2.53 ± 1.21 


H 14- 


-^9 


12.587 


12.588 ± 0.0024 


6.90 ± 1.996 




2.91 ± 1.11 


H 11- 


-^8 


12.387 




Blended with 7^6 




H7^ 


6 


12.372 


12.364 ± 0.0040 


12.49 ± 3.691 




4.41 ± 1.71 


H20- 


-^10 


12.157 


12.155 ± 0.0048 


12.20 ± 4.486 




2.88 ± 1.38 


H15- 


-^9 


11.540 


11.536 ± 0.0020 


7.88 ± 1.774 




3.29 ± 0.97 


H9^ 


7 


11.309 


11.305 ± 0.0024 


69.90 ± 2.115 




2.38 ± 0.95 


H23- 


-^10 


11.243 


11.238 ± 0.0039 


6.00 ± 3.846 




1.33 ± 1.04 


H24- 


-^10 


11.033 


11.033 ± 0.0052 


9.45 ± 4.707 




2.20 ± 1.47 


H 16- 


-^9 


10.804 


10.801 ± 0.0019 


7.68 ± 1.812 




4.24 ± 1.35 


H26- 


-^10 


10.701 


10.698 ± 0.0036 


5.10 ± 2.802 




1.63 ± 1.17 


H27- 


-^10 


10.567 


10.577 ± 0.0058 


6.70 ± 4.402 




1.17 ± 1.05 


H 12- 


-^8 


10.504 


10.502 ± 0.0023 


6.04 ± 2.163 




2.55 ± 1.27 


H28- 


-^10 


10.451 


10.452 ± 0.0067 


4.690 ± 6.250 




7.24 ± 1.22 


H 17- 


-^9 


10.261 


10.258 ± 0.0018 


7.51 ± 1.642 




4.47 ± 1.29 


H30- 


-^10 


10.258 




Blended with 17^9 




H 18- 


-^9 


9.847 


9.845 ± 0.0017 


5.29 ± 1.330 




5.25 ± 1.56 


H 19- 


-^9 


9.522 


9.518 ± 0.0038 


5.59 ± 3.424 




2.92 ± 2.07 


H 13- 


-^8 


9.392 


9.391 ± 0.0017 


5.27 ± 1.597 




4.55 ± 1.79 


H20- 


-^9 


9.261 


9.257 ± 0.0026 


5.67 ± 2.591 




3.49 ± 2.02 


H21- 


-^9 


9.047 


9.047 ± 0.0022 


5.20 ± 2.006 




3.49 ± 1.76 


H22- 


-^9 


8.870 


8.871 ± 0.0073 


10.41 ± 7.766 




6.53 ± 5.99 


H 10- 


-^7 


8.760 


8.758 ± 0.0043 


5.43 ± 5.002 




3.40 ± 4.26 


H29- 


-^9 


8.173 




Blended with 15^8 




H 15- 


-^8 


8.155 


8.152 ± 0.0014 


3.81±1.021 




7.0299±2.5343 
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Table 2. Data on observed lines 









7 ' 


Cas, Aug 1998 




Transition 


■^theory 


^obs 


AA 


Fu 


H7^ 


6 


12.372 


12.369 ± 0.0031 


9.53 ± 2.643 


1.09 ± 0.363 


H9^ 


7 


11.309 


11.307 ± 0.0015 


6.57 ± 1.391 


1.20 ± 0.328 


H 10- 


-^7 


8.760 


8.759 ± 0.0031 


6.13 ± 2.804 


2.19 ± 1.328 


H 11- 


-^8 


12.387 


12.394 ± 0.0065 


7.71 ± 6.266 


0.45 ± 0.433 


H 12- 


-^8 


10.504 


10.504 ± 0.0018 


6.85 ± 1.670 


1.25 ± 0.400 


H 13- 


-^8 


9.392 


9.393 ± 0.0015 


6.85 ± 1.386 


1.79 ± 0.474 


H 14- 


-^8 


8.665 


8.666 ± 0.0028 


7.69 ± 2.543 


2.33 ± 1.174 


H 14- 


-^9 


12.587 


12.589 ± 0.0022 


6.37 ± 1.716 


1.05 ± 0.370 


H 15- 


-^8 


8.155 


8.154 ± 0.0010 


4.20 ± 0.832 


4.96 ± 1.311 


H 15- 


-^9 


11.540 


11.540 ± 0.0017 


8.03 ± 1.555 


1.11 ± 0.283 


H 16- 


-^9 


10.804 


10.803 ± 0.0015 


6.39 ± 1.314 


1.41 ± 0.383 


H 17- 


-^9 


10.261 


10.260 ± 0.0018 


7.62 ± 1.620 


1.28 ± 0.361 


H 18- 


-^9 


9.847 


9.846 ± 0.0016 


6.38 ± 1.349 


1.74 ± 0.494 


H 18- 


-^10 


13.188 


13.197 ± 0.0034 


6.70 ± 2.798 


0.74 ± 0.410 


H 19- 


-^9 


9.522 


9.522 ± 0.0023 


2.70 ± 2.071 


0.74 ± 0.721 


H 19- 


-^10 


12.611 


12.615 ± 0.0021 


5.63 ± 2.998 


0.87 ± 0.615 


H20- 


-^9 


9.261 


9.262 ± 0.0035 


9.17 ± 3.226 


0.88 ± 0.433 


H20- 


-^10 


12.157 


12.157 ± 0.0044 


9.53 ± 4.018 


0.46 ± 0.259 


H21- 


-^9 


9.047 


9.048 ± 0.0017 


2.27 ± 2.738 


2.17 ± 5.413 


H21- 


-^10 


11.792 


11.793 ± 0.0043 


7.66 ± 3.738 


0.52 ± 0.334 


H22- 


-^9 


8.870 


8.873 ± 0.0057 


6.64 ± 5.096 


1.34 ± 1.340 


H22- 


-^10 


11.492 


11.488 ± 0.0032 


4.14 ± 4.405 


0.36 ± 0.537 


H23- 


-^9 


8.721 


8.726 ± 0.0172 






H23- 


-^10 


11.243 


11.241 ± 0.0054 


7.70 ± 4.900 


0.39 ± 0.334 


H24- 


-^9 


8.594 


8.594 ± 0.0021 


2.10 ± 2.184 


2.02 ± 3.400 


H24- 


-^10 


11.033 


11.034 ± 0.0061 


8.01 ± 5.466 


0.43 ± 0.383 


H25- 


-^9 


8.485 


8.489 ± 0.0050 


6.74 ± 4.387 


1.43 ± 1.237 


H25- 


-^10 


10.855 


10.849 ± 0.0099 






H26- 


-^9 


8.391 


8.390 ± 0.0036 


4.28 ± 3.611 


1.61 ± 1.737 


H26- 


-^10 


10.701 


10.694 ± 0.0096 






H27- 


-^9 


8.309 


8.307 ± 0.0053 






H27- 


-^10 


10.567 


10.580 ± 0.0091 






H28- 


-^9 


8.236 


8.240 ± 0.0018 


3.65 ± 1.244 


3.58 ± 1.550 


H28- 


-^10 


10.451 


10.453 ± 0.0248 






H29- 


-^9 


8.173 


Blended with 15^8 


H29- 


-^10 


10.348 


10.334 ± 0.0049 






H30- 


-^9 


8.116 


8.107 ± 0.0040 


6.80 ± 3.654 


1.71 ± 1.218 


H30- 


-^10 


10.258 


Blended with 17^9 


H31- 


-^10 


10.177 


10.181 ± 0.0074 


6.71 ± 6.115 


0.28 ± 0.332 
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Table 2G. Data on 


observed lines - 


o Aqr 




Transition 


■^theory 


Ao6s 


AA 




F, 






(^m) 


(^m) 


(10"-Vm) 


(10- 


"^^ W cm ^ ^m. ^) 


H 19- 


-^10 


12.611 


12.618 ± 0.0162 


9.24 ± 13.21 




2.03 ± 19.08 


H 14- 


-^9 


12.587 


12.590 ± 0.0041 


10.84 ± 3.713 




2.78 ± 1.23 


H 11- 


-^8 


12.387 


12.384 ± 0.0084 


21.00 ± 7.290 




4.24 ± 1.98 


H7^ 


6 


12.372 


12.370 ± 0.0027 


5.44 ± 1.812 




1.86 ± 0.83 


H21- 


-^10 


11.792 


11.782 ± 0.0061 


12.13 ± 5.769 




1.43 ± 0.94 


H15- 


-^9 


11.540 


11.543 ± 0.0035 


9.75 ± 3.229 




18.11 ± 6.02 


H9^ 


7 


11.309 


11.307 ± 0.0020 


9.35 ± 1.813 




2.68 ± 0.69 


H24- 


-^10 


11.033 


11.024 ± 0.0081 








H 16- 


-^9 


10.804 


10.805 ± 0.0036 


7.73 ± 3.211 




1.33 ± 0.72 


H 12- 


-^8 


10.504 


10.502 ± 0.0013 


7.26 ± 1.142 




2.89 ± 0.60 


H28- 


-^10 


10.451 


10.453 ± 0.0264 








H29- 


-^10 


10.348 


10.343 ± 0.0020 








H 17- 


-^9 


10.261 


10.260 ± 0.0015 


4.80 ± 1.005 




1.72 ± 0.48 


H30- 


-^10 


10.258 




Blended with 17^9 




H 18- 


-^9 


9.847 


9.846 ± 0.0029 


4.70 ± 2.315 




1.33 ± 0.81 


H 13- 


-^8 


9.392 


9.392 ± 0.0023 


8.80 ± 2.075 




3.30 ± 1.04 


H22- 


-^9 


8.870 


8.871 ± 0.0116 








H 10- 


-^7 


8.760 


8.759 ± 0.0026 


7.86 ± 2.421 




4.61 ± 1.86 


H 14- 


-^8 


8.665 


8.668 ± 0.0024 


5.10 ± 2.145 




3.00 ± 1.67 


H 15- 


-^8 


8.155 


8.155 ± 0.0011 


4.43 ± 1.011 




4.68 ± 1.40 


H29- 


-^9 


8.173 




Blended with 15^8 
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Table 2H. Data on observed lines - EW Lac 



Transition 


■^theory 


^obs 


AA 




Fu 






(/wm) 


(//m) 


(10-3/im) 


(10- 


-19 W cm-2 ^m^i) 


H 19- 


-^10 


12.611 


12.614 ± 0.0076 


8.70 ± 6.928 




1.24 ± 1.16 


H 14- 


-^9 


12.587 


12.591 ± 0.0045 


8.40 ± 3.521 




2.03 ± 1.05 


H 11- 


-^8 


12.387 




Blended with H 7^' 


6 


H7^ 


6 


12.372 


12.373 ± 0.0017 


12.40 ± 1.549 




3.73 ± 0.62 


H15- 


-^9 


11.540 


11.543 ± 0.0026 


8.20 ± 2.307 




0.88 ± 0.70 


H22- 


-^10 


11.492 


11.491 ± 0.0100 








H9^ 


7 


11.309 


11.308 ± 0.0017 


8.40 ± 1.575 




1.66 ± 0.41 


H 16- 


-^9 


10.804 


10.805 ± 0.0023 


6.60 ± 1.876 




0.68 ± 0.27 


H26- 


-^10 


10.701 


10.696 ± 0.0029 


4.83 ± 3.181 




0.39 ± 0.32 


H 12- 


-^8 


10.504 


10.500 ± 0.0021 


7.60 ± 1.871 




1.09 ± 0.36 


H 17- 


-^9 


10.261 


10.263 ± 0.0027 


14.00 ± 2.669 




1.78 ± 1.16 


H30- 


-^10 


10.258 




Blended with 17^S 


1 


H 19- 


-^9 


9.522 


9.523 ± 0.0024 


5.50 ± 2.010 




1.41 ± 1.44 


H 13- 


-^8 


9.392 


9.386 ± 0.0026 


5.70 ± 2.528 




1.11 ± 1.39 


H20- 


-^9 


9.261 


9.254 ± 0.0043 


3.86 ± 2.720 




5.50 ± 4.00 


H21- 


-^9 


9.047 


9.048 ± 0.0027 


2.50 ± 3.141 




3.96 ± 5.02 


H 10- 


-^7 


8.760 


8.758 ± 0.0014 


6.10 ± 1.321 




1.79 ± 1.11 


H 14- 


-^8 


8.665 


8.664 ± 0.0031 


6.60 ± 2.555 




0.85 ± 0.97 
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Table 21. Data on observed lines - n Dra 



Transition 


■^theory 


^obs 


AA 


Fu 






(^m) 


(fim) 


(10-3;^m) (10- 


19 w cm-2 /xm'i) 


H 18- 


-^10 


13.188 


13.192 ± 0.0069 


6.30 ± 7.434 


1.67 ± 2.47 


H 19- 


-^10 


12.611 


12.617 ± 0.0026 


7.50 ± 2.327 


1.34 ± 1.73 


H 14- 


-^9 


12.587 


12.590 ± 0.0041 


8.40 ± 3.497 


2.56 ± 2.30 


H 11- 


-^8 


12.387 


12.393 ± 0.0026 


3.54 ± 3.874 


0.18 ± 2.83 


H7^ 


6 


12.372 


12.376 ± 0.0017 


7.55 ± 1.727 


0.35 ± 1.75 


H20- 


-^10 


12.157 


12.155 ± 0.0046 


6.93 ± 3.672 


1.07 ± 1.73 


H21- 


-^10 


11.792 


11.793 ± 0.0059 


12.70 ± 5.418 


1.96 ± 2.42 


H15- 


-^9 


11.540 


11.539 ± 0.0021 


10.30 ± 1.897 


3.03 ± 1.61 


H22- 


-^10 


11.492 


11.491 ± 0.0055 


8.00 ± 4.943 


0.90 ± 0.75 


H9^ 


7 


11.309 


11.309 ± 0.0014 


7.60 ± 1.341 


2.75 ± 1.35 


H23- 


-^10 


11.243 


11.244 ± 0.0061 


7.80 ± 5.390 


0.87 ± 1.76 


H25- 


-^10 


10.855 


10.847 ± 0.0042 


8.01 ± 12.14 


1.16 ± 2.33 


H 16- 


-^9 


10.804 


10.803 ± 0.0014 


79.40 ± 1.301 


3.39 ± 1.61 


H26- 


-^10 


10.701 


10.701 ± 0.0027 






H 12- 


-^8 


10.504 


10.504 ± 0.0017 


8.00 ± 1.520 


3.10 ± 1.70 


H 17- 


-^9 


10.261 


10.263 ± 0.0017 


8.50 ± 1.528 


3.33 ± 1.74 


H 18- 


-^9 


9.847 


9.845 ± 0.0022 


6.92 ± 1.979 


2.71 ± 1.03 


H 19- 


-^9 


9.522 


9.521 ± 0.0023 


8.57 ± 2.067 


3.49 ± 1.11 


H 13- 


-^8 


9.392 


9.390 ± 0.0014 


6.75 ± 1.202 


3.84 ± 0.93 


H20- 


-^9 


9.261 


9.260 ± 0.0034 


5.20 ± 2.654 


12.16 ± 6.23 


H21- 


-^9 


9.047 


9.047 ± 0.0033 


7.82 ± 2.956 


2.59 ± 1.27 


H 10- 


-^7 


8.760 


8.760 ± 0.0012 


5.75 ± 1.054 


3.81 ± 0.93 


H 14- 


-^8 


8.665 


8.666 ± 0.0018 


7.48 ± 1.604 


3.92 ± 1.12 


H24- 


-^9 


8.594 


8.592 ± 0.0026 






H15- 


-^8 


8.155 


8.158 ± 0.0009 


6.49 ± 0.814 


6.58 ± 1.09 


H29- 


-^9 


8.173 




Blended with 15^8 
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Table 3A. Data on observed lines - (3 Lyr 



Transition^ 


^theory 


^obs 


AA 


F, 




(^m) 


(^m) 


(10-Vm) 


(10-19 W cm-2 ^m-i) 


[Ne II] 


12.81 


12.813 ± 0.0017 


9.91 ± 1.565 


8.83 ± 1.84 


H 14^9 


12.587 


12.583 ± 0.0039 


10.60 ± 3.615 


5.08 ± 2.26 


H 11^8 


12.387 


12.393 ± 0.0035 


6.84 ± 3.384 


4.09 ± 2.31 


H7^6 


12.372 


12.368 ± 0.0020 


10.16 ± 1.851 


38.90 ± 8.70 


H9^7 


11.309 


11.305 ± 0.0004 


9.86 ± 0.389 


20.81 ± 1.09 


He I 4P-4S, 3PO-3S 


10.882 


10.882 ± 0.0010 


8.69 ± 0.953 


14.24 ± 2.07 


H 16^9 


10.804 


10.801 ± 0.0071 


6.85 ± 6.682 


2.18 ± 2.65 


H 12^8 


10.504 


10.504 ± 0.0014 


10.55 ± 1.306 


14.23 ± 2.34 


H 17^9 


10.261 


10.263 ± 0.0066 


7.37 ± 5.992 


2.45 ± 2.58 


H 18^9 


9.847 


9.846 ± 0.0028 


8.62 ± 2.576 


3.12 ± 1.25 


? 




9.597 ± 0.0012 


4.34 ± 1.562 


2.49 ± 1.12 


H 19^9 


9.522 


9.516 ± 0.0027 


5.74 ± 2.232 


2.70 ± 1.44 


H 13^8 


9.392 


9.389 ± 0.0012 


9.77 ± 1.029 


10.31 ± 1.43 


H 10^7 


8.760 


8.758 ± 0.0007 


8.14 ± 0.608 


23.92 ± 2.36 


H 14^8 


8.665 


8.663 ± 0.0030 


7.64 ± 2.760 


7.68 ± 3.53 


H 15^8 


8.155 


8.153 ± 0.0020 


5.11 ± 1.780 


6.37 ± 2.75 



"^The hydrogen recombination lines may also be blended with the identical He I recombi- 
nation transitions. 
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Table 3B. Data on observed lines - MWC 349 



Transition Xtheory Kbs AA F,, 

(fiui) (^m) (10~^/im) (lO^-*^^ W cm~^ /im" 



[Ne II] 


12.81 


12.813 ± 0.0002 


9.55 ± 0.168 


95.35 ± 2.21 


H7^6 


12.372 


12.372 ± 0.0017 


13.35 ± 1.534 


10.80 ± 1.64 


? 




11.918 ± 0.0045 


14.39 ± 4.390 


5.89 ± 2.32 


H9^7 


11.309 


11.310 ± 0.0008 


7.83 ± 0.771 


0.82 ± 0.70 


H 12^8 


10.504 


10.508 ± 0.0025 


9.98 ± 2.329 


0.06 ± 1.18 
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Table 4: Properties of our Sources 



Object 


Shell 


Spectroscopic 


X-Ray 


Radio 


References 




Star 


Binary 


Binary 


Source 




ip Per 


yes 


no 


no 


yes 


1,7 


rj Tau 


yes 


yes 


no 


yes 


1,3,9 


CTau 


yes 


yes 


no 


no 


1, 3 


48 Per 


no 


yes 


yes 


no 


2,6 


(/)Per 


yes 


yes 


no 


no 


1 


7 Cas 


yes 


no 


yes 


yes 


1,4,9 


o Aqr 


yes 


no 


no 


no 


1 


EWLac 


yes 


no 


no 


yes 


1,9 


K Dra 


no 


yes 


yes 


no 


5, 6 


/3Lyr 


no 


yes 


no 


yes 


5, 10 


MWC 349 


no 


no 


no 


yes 


8 



(1) Slettebak et al. 1992, (2) Poeckert 1981, (3) Jared et al. 1989, (4) Haberl 1995, (5) Pols et al. 1991, (6) Peters 
1982, (7) Marlborough et al. 1997, (8) Strelnitski et al. 1996, (9) Taylor et al. 1990, (lO)Pesek et al 1996. 
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Table 5. Stellar Parameters for Sources 



Object 


^m 


R* 


D 


logL 


1.25Atm 


2.2^m 




(Ro) 


(pc) 


(logLo) 


(mag) 


(mag) 


V' Per 


16000 


7.0 


215 


3.20 


4.30 


3.94 


rj Tau 


14000 


6.5 


113 


2.86 


2.95 


2.96 


C Tau 


16500 


6.5 


128 


3.05 


3.20 


2.94 


48 Per 


16000 


6.3 


170 


3.60 


3.99^ 


3.79^ 


(/^Per 


20000 


7.5 


220 


3.20 


3.83 


3.25 


7 Cas 


31000 


10.0 


188 


4.93 


2.47 


2.20 


o Aqr 


11000 


4.0 


117 


2.45 


4.84 


4.80 


EWLac 


16000 


6.5 


337 


3.54 


5.34 


5.00 


K Dra 


14500 


6.4 


153 


2.70 


3.86 


3.75 



'"At l.l;um. 
^At 2.3^m. 
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Table 6. Parameter Values for Free-Free Fits 



Object 




Tsh 




rif. 


Rsh 


X^ 




(] 


LO^K) 


(10^ 


1 cm-3) 


(R.) 




■0 Per 


1.60 


;+.07/-.15) 


2.60 


+.16/-.07) 


4.61 (+.21/-.20) 


1.56 


ryTau 


0.80 


'+.02/-.06) 


1.52 


+.08/-.03) 


3.07 (±.15) 


1.04 


CTau 


1.39 


'+.04/-.10) 


2.30 


+.11/-.04) 


4.95 (±.32) 


4.27 


48 Per 


1.29 


'+.05/-.08) 


2.21 


+.09/-.05) 


4.04 (±.13) 


1.05 


(pFei 


1.65 


'+.06/-.08) 


3.01 


+.19/-.10) 


5.02 (±.10/-.ll) 


0.74 


7 Cas 


2.68 


'+.15/-.14) 


4.57 


+.33/-.21) 


6.00 (±.13) 


1.12 


Aqr 


1.01 


:+.02/-.ll) 


1.60 


+.05/-.03) 


2.99 (±.11/-.13) 


0.54 


EWLac 


1.32 


'+.02/-.09) 


2.20 


+.08/-.03) 


3.80 (±.20/-.22) 


1.49 


K Dra 


1.20 


:+.o4/-.ii) 


1.83 


+.14/-.04) 


3.12 (±.14/-.16) 


2.13 
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Table 7. Comparision of Free-Free Fits 



Object 


Author"" 


n 


R. 


Tsh 


Rsh 


rie 


Mh 






(10^ K) 


(Rq) 


(10^ K) 


(Ro) 


(10^^ cm-3) 


(10-10 Mo) 


ip Per 


RHS 


1.6 


7.0 


1.6 


32.3 


2.60 


9.3 




WCL 


L5 


6.0 


1.2 


37.2 






7?Tau 


RHS 


1.4 


6.5 


0.8 


20.0 


1.52 


0.93 


C Tau 


RHS 


1.65 


6.5 


1.39 


32.2 


2.30 


8.3 




GHJ 


1.8 


5.75 


1.4 


20.1 


4.4 






WCL 


1.4 


6.0 


1.1 


28.8 






48 Per 


RHS 


1.6 


6.3 


1.29 


25.5 


2.21 


3.7 




WCL 


1.8 


6.5 


1.45 


36.4 






(j)FeT 


RHS 


2.0 


7.5 


1.65 


37.7 


3.01 


17.6 




GHJ 


1.6 


5.0 


1.2 


23.0 


3.8 






WCL 


2.0 


7.0 


1.6 


91.0 






7 Cas 


RHS 


3.1 


9.0 


2.68 


54.0 


4.57 


81.4 




GHJ 


2.6 


7.2 


2.0 


29.5 


4.4 






WCL 


3.1 


10.0 


2.5 


80.0 






o Aqr 


RHS 


1.1 


4.0 


1.01 


12.0 


1.9 


0.25 




WCL 


1.3 


3.4 


1.05 


9.5 






EWLac 


RHS 


1.6 


6.5 


1.32 


24.7 


2.20 


3.2 




WCL 


1.4 


5.7 


1.1 


27.4 






K Dra 


RHS 


1.45 


6.5 


1.20 


20.3 


1.83 


1.2 




WCL 


1.3 


4.5 


1.05 


22.1 







^WCL = Waters et al. 1987, GHJ = Gehrz et al 1974, and RHS = this work. 
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Fig. 1. — The reduced spectra of our sources, (a) ^ Per, (b) rj Tau, (c) (" Tau, (d) 48 Per, (e) (p Per, 
(f) 7 Cas. 
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Fig. 1, Continued. — The reduced spectra of our sources, continued, (g) o Aqr, (h) EW Lac, (i) 
K Dra, (j) (3 Lyr, and (k) MWC 349. 
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Fig. 2. — The spectrum produced by coadding all nine of the hydrogen spectra (HS) Be star spectra. 
Several lines which are visible in the individual spectra become clearer in this summed spectra. 
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Fig. 3. — A comparison plot of the H 9^7/H 7^6 ratio versus the H 10^7/H 7— >6 ratio. No 
strong correlation can be seen in this small sample set. 
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Fig. 4. — The spectrum of 7 Cas (dotted line), with the free- free emission fit overlaid (solid line). 
The dashed line and dash-dot line are the optically thin and optically thick components of the 
emission, respectively. The open diamonds are the points used to make the fit, and the open 
triangle is the IRAS 12/im flux. 



